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Abstract: In this talk we present a calculation of the dacay rate difference in the neutral Bg — Bg 
system, AF^^, in next-to- leading order (NLO) QCD. We find a sizeable decrease compared to lead- 
ing-order (LO) estimates: (Ar/r)^, = (/i3,/210Mel/)2[0.006B(mb)-|-0.I50Bs(TO6)-0.0063] in terms 
of the bag parameters B and Bs in the NDR scheme. We put special emphasize on the theoretical 
and physical implications of this quantity. 



1. Non-expert-introduction 

As there were many students in the audience we 
will start with an elementary introduction. Neu- 
tral mesons are well known from lectures at the 
university and were mentioned here several times 
e.g. in 0, g, ^, Q|. As in the K-system we have 
in the Ss-system flavour eigenstates which are 
defined by their quark content. 



\Bg) = (bs) ; \Bs) = (bs) . 



(1.1) 



The mass eigenstates are linear combinations of 
the flavour eigenstates 



\Bh) = p\Bs 
\Bl) - p\Bs 



q\Bs 
q\Bs 



(1.2) 
(1.3) 



1. 



with the normalization condition |pp -I- \q 
Bh and Bl are the physical states. They have 
definite masses and lifetimes, but no definte CP- 
quantum numbers. The mass eigenstates are in 
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general mixtures of CP-odd and CP-even eigen- 
states. 

The time evolution of the physical states is 
described by a simple Schrodinger equation 



idtB = HB 



(1.4) 



with 



LB,, 



H = 



M-^ 



Mt 



2^ 12 



Mi2 - 

M ■ 



5^12 



(1.5) 

To find the mass eigenstates and the eigenvalues 
of the mass operator and the dacay rate operator 
we have to diagonalize the hamiltonian. We get 



AAfs - Mh - Ml = 2Re{Q) 
ATb - Fl - Fh = 4/m(0) 



with 



0-\/(M: 



-ri2)(M: 



12 



12; 



(1.6) 
(1.7) 



(1.8) 



If we neglect CP violation and expand in ra^/ra^ 
we can write with a very good precision 



AAfs = 2IM12I 
AFb = -2Fi2 



(1.9) 
(1.10) 



Corfu Summer Institute on Elementary Particle Physics, 1998 



Alexander Lenz 



The different neutral meson systems gave rise to 
important contributions to the field of high en- 
ergy physics. In 1964 Christenson, Cronin, Fitch 
and Turlay [p| discovered indirect CP-violation^ 
in the K^ — K -system. The mass difference 
in the Bd — Bd -system was the first experimen- 
tal hint for a very large top quark mass, before 
the indirect determination at LEP and before the 
discovery at Tevatron. As mt is by now quite 
well known, we can extract the CKM parameter 
|Vtd^h| from AMb^- The determination of the 
CKM parameters is crucial for a test of our un- 
derstanding of the standard model and for the 
search for new physics. The mass difference in 
the Bs — Bs -system is not measured yet, but we 
have a lower limit from which we already get an 
important bound on the parameters of the CKM 
matrix. 

The Heavy Quark Expansion (HQE) is the the- 
oretical framework to handle inclusive B-decays. 
It allows us to expand the dacay rate in the fol- 
lowing way 



mfc / \nn, 



(1.11) 



Here we have an systematic expansion in the 
small parameter A /mi,. The different terms have 
the following physical interpretations: 

• To'. The leading term is described by the 
decay of a free quark (parton model), we 
have no non-perturbative corrections. 



• Fi: In the derivation of eq. (1.11) we make 
an operator product expansion. From di- 
mensional reasons we do not get an opera- 
tor which would contribute to this order in 
the HQE. 2 

• F2: First non-perturbative corrections arise 
at the second order in the expansion due to 
the kinetic and the chromomagnetic opera- 
tor. They can be regarded as the first terms 
in a non-relativistic expansion. 



^Indirect or equivalently mixing induced CP violation 
means that the physical states Kg/j^ are not pure CP- 
eigenstates. There is a big contribution of one CP-parity 
and a tiny of the opposite CP-parity. If the small contri- 
bution dacays, one speaks of indirect CP violation. 

^Strictly spoken we get one operator of the appropri- 
ate dimension, but with the equations of motion we can 
incorporate it in the leading term. 



• Fa: In the third order we get the so-called 
weak annihilation and pauli interference di- 
agrams. Here the spectator quark is in- 
cluded for the first time . These diagrams 
give rise to lifetime differences in the neu- 
tral i3-system. 

Each of these terms can be expanded in a power 
series in the strong coupling constant 



r, ^ Fi°) + ^FP) 



So ATb has the following form 



(1.12) 



Ar^ . 4 (rf + ^F« 

mt \ TT 



A' 



(rf + ;jr.-. + ..>^(r<". + ...) 

(1.13) 

After this short introduction for non-experts we 
motivate the special interest in the quantity AF^, . 



2. Motivation 

From a physical point of view one wants to know 
the exact value of the decay rate difference, be- 
cause 

• (AF/F)b^ is expected to be large. LO §] 
estimates give values up to 20%. This is 
on the border of the experimental visibility 

0; 

• a big value of AF^^ would enable us to do 
novel studies of CP-violation without the 
need of tagging |^. Tagging is a major ex- 
permintal difficulty in B-physics; 

• in the ratio AF^^ /AMb^ some of the non- 
perturbative parameters cancel 



10. So 



we can get theoretically clean information 
on AMb^ from a measurement of AF^^ ; 

• the decay rate difference can be used to 
search for non SM-physics. In |11 it was 

shown that AT new physics < ATsM. 

In order to fuUfiU this physics program we need a 
relieable prediction in the standard model. There- 



Co) 
3 ' 



fore we need in addition to the LO estimate F 
which was calculated in [p| 

• the 1/mb-corrections F4 . They have been 
calculated by M; 
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• the non-perturbative matrix elements for 
the Ai? = 2 operators, which arise in the 
calculation. Here a relieable prediction is 
still missing; 

• the NLO QCD corrections to the leading 
term in the 1/mi, expansion, Fg . This was 
the aim of our work |lO| . Besides the bet- 
ter accuracy and a reduction of the /i de- 
pendence there is a very important point: 
NLO-QCD correction are needed for the 
proper matching of the perturbative calcu- 
lation to lattice calculations. 

From a technical point of view this calculation 
was very interesting because 

• our result provides the first calculation of 
perturbative QCD corrections beyond lead- 
ing logarithmic order to spectator effects 
in the HQE. Soft gluon emmision from the 
spectator s quark leads to power-like in- 
frared singularities in individual contribu- 
tions. As a conceptual test of the HQE the 
final result has to be infrared finite |lj] . 

• a crucial point in the derivation of the HQE 
is the validity of the operator product ex- 
pansion. This assumption is known under 
the name quark hadron duality and can be 
tested via a comparison of theory and ex- 
periment. A recent discussion of that sub- 
ject can be found in Q. 

In the next chapter we will describe the calcula- 
tion. 

3. Calculation 

The width difference in the B^ — B -system is 
defined as 



AF = Fr -F 



H 



-2F21 



(3.1) 



The off-diagonal element of decay-width matrix 
can be related to the so-called transition operator 
T via 

r2i = 77^{Bs\T\Bs) (3.2) 



2M 



Bs 



with 



T = Imi I d^xT Heff {x)Heff (0) . (3.3) 



In T we have a double insertion of the effective 
hamiltonian with the standard form ||14| 



n 



eff 



^V:,V,JY,CrQr + CsQs] . (3.4) 



Gf denotes the Fermi constant, Vpq are the CKM 
matrix elements and Qi are local AB = 1 oper- 
ators. The Wilson coefficients Ci describe the 
short distance physics and are known to NLO 
QCD. 

Formally we proceed now with an operator prod- 
uct expansion of that product of two hamiltoni- 
ans. In real life one has to calculate diagrams of 
the following form: 




Figure 1: The imaginary part of massive two loop 
diagram of that form has to be calculated. 

One can do the calculation in two different ways 
(we did it in both ways, to have a check): 

• calculate the imaginary part of the two loop 
integrals 



• use Cutkosky rules and calculate virtual 
and real one loop corrections, followed by 
a phase space integration. 

The result in LO QCD has the following form 

^2 ^2 



r 



-%^ iV;,V,sf [G{z)Q + Gs{z)Qs] 



(3.5) 



with z = rri^/ml and the /S.B — 2 operators 



Q = {biSi)v-A{hjSj)v-A 
Qs = (biSi)s-pibjSj)s^p ■ 



(3.6) 



In principle we have more operators, but we can 
reduce them to the two operators above with the 
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use of Fierz identities 



Equation ( 3^ ) is an example of an operator prod- 
uct expansion of equation (3_^). We have re- 
duced the double insertion of Ai? = 1 operators, 
which appear in tteff, to a single insertion of 
an AB = 2 operator. In principle we have inte- 
grated out the internal charm quarks in figure 1. 
For the NLO calculation we have to match the 
A_B = 1 double insertion with gluon exchange to 
a /S.B = 2 insertion with gluon exchange. This 
means, we have to calculate the following dia- 
grams: 




\o^ '^^^5t 



Figure 2: All diagrams that have to be calculated 
for the NLO QCD determination of AVb, ■ 

These diagrams can be classified in the fol- 
lowing way: 

El — E^: Virtual one loop corrections to a A_B = 2 
operator insertion. 

Di — DiQ : Imaginary part of virtual two loop correc- 
tions to a double insertion of AB = 1 op- 
erators. 

-Dii, D12: Penguin contributions to the AB — 1 dou- 
ble insertion. 

The calculation of all these diagrams gives us the 
NLO QCD result. 

4. Results 

The result in NLO is: 



T 



Gj^rUh ,,,»,, ,2 



127T 



{V;,V,,r[G{z)Q-Gs{z)Qi 



(4.1) 



'This reduction is relativcley tricky. For details see 



with the following numerical values for the Wil- 
son coefficients 



Ai 


mfc/2 


TOfc 


2mb 


G(0) 


0.013 


0.047 


0.097 


G 


0.023 


0.030 


0.036 


gT 


1.622 


1.440 


1.292 


Gs 


0.743 


0.937 


1.018 



with 



G = G(") 



47r 



(4.2) 



Here one can see two important points. First, the 
value for Gs is numerical dominant and second, 
the NLO values are considerably smaller than the 
LO values. 

For the final result we parametrise the ma- 
trix elements of the AB — 2 operators in the fol- 
lowing way: 



{Bs\Q\B, 
{B,\Qs\B, 



-jlMlB 



Ml 



3^^=^^^=(m, + m.)2 



Bs 



B and Bs are so-called bag parameters, f^^ is 
the decay constant. The values of these param- 
eters have to be determined by non-perturbative 
methods like lattice simulations, rhq denotes the 
running quark mass in the 7\f 5-scheme. 
With the following input parameters 



m,b = 4.8GeV 



nib 



0.085 fhs = 0.2GeV 



Ms, = 5.37GeV B{Bs 1-^ Xev) = 0.104 

(4.3) 
we obtain for the relative dacay rate difference 



( r )b, - VsioMeVy 



fE 



V 



[0.006B{mb) + 0.l50Bsimb) - 0.063] 



(4.4) 
A definitive determination of the two bag param- 
eters is still missing. From the literature p5| we 
were able to extract preliminary values for the 
bag parameters 



B{mb) = 0.9 Bsirrib) = 0.75 



(4.5) 
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With that numbers at hand we obtain as a final 
result 



References 



Ar 



/, 



Bs 



210MeV 



(0.054ig;°i^±???) 



1.032- 



(4.6) 

The question marks remind us that we do not 
know the uncertainties in the numerical values 
for the bag parameters. 

5. Disscussion and outlook 

The LO estimate for the relative decay rate dif- 
ference AFs^/Tb^ = 0(20%) is considerably re- 
duced due to several effects: 

• the l/rrifc corrections are sizeable and give 
an absolute reduction of about - 6.3 % Q. 

• the pure NLO QCD corrections are size- 
able, too and give an absolute reduction of 
about - 4.8 % 0- 

• with the NLO QCD corrections at hand we 
can perform a proper matching to the (pre- 
liminary) lattice calculations for the bag 
parameters. This tells us that we have to 
use a low value for the bag parameters, i.e. 
Bsinn,) = 0.75 0, [5|. Compared to the 
naive estimate Bs — 1, this is another ab- 
solute reduction of about - 3.8 %. 

Unfortunateley the value of AF^^ /^b^ has been 
pinned down to a value of about 5%. The LO 
prediction was just a the border of experimental 
visibility |Q. Now we will have to wait for the 
forthcoming experiments like HERA-B, Tevatron 
(run II) and LHC. 

Another application of our calculation are in- 
clusive indirect CP-asymmetries in the b — > uud 
channel. For the complete NLO prediction of 
this quantity, F12 in the Bg system was missing. 



We get this value from our calculation with a 



trivial exchange of the CKM parameters and the 
limit nic 1— *■ 0. This allows a determination of the 
CKM-anglea 0,|l|. 
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